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OBJECTIVES: This study evaluated plantar thermography sensitivity and specificity in diagnosing diabetic 
polyneuropathy using cardiac tests (heart rate variability) as a reference standard because autonomic small fibers 
are affected first by this disease. 

METHODS: Seventy-nine individuals between the ages of 19 and 79 years old (28 males) were evaluated and divided 
into three groups: control (n = 37), pre-diabetics (n = 13) and type 2 diabetics (n = 29). The plantar images were 
recorded at baseline and then minutes after a provocative maneuver (Cold Stress Test) using an infrared camera that 
is appropriate for clinical use. Two thermographic variables were studied: the thermal recovery index and the 
interdigital anisothermal technique. Heart rate variability was measured in a seven-test battery that included three 
spectral indexes (in the frequency domain) and four Ewing tests (the Valsalva maneuver, the orthostatic test, a deep 
breathing test, and the orthostatic hypotension test). Other classically recommended tests were applied, including 
electromyography (EMG), Michigan inventory, and a clinical interview that included a neurological physical 
examination. 

RESULTS: Among the diabetic patients, the interdigital anisothermal technique alone performed better than the 
thermal recovery index alone, with a better sensitivity (81.3%) and specificity (46.2%). For the pre-diabetic patients, 
the three tests performed equally well. None of the control subjects displayed abnormal interdigital anisothermal 
readouts or thermal recovery indices, which precluded the sensitivity estimation in this sample of subjects. However, 
the specificity (70.6%) was higher in this group. 

CONCLUSION: In this study, plantar thermography, which predominately considers the small and autonomic fibers 
that are commonly associated with a sub-clinical condition, proved useful in diagnosing diabetic neuropathy early. 
The interdigital anisothermal test, when used alone, performed best. 
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INTRODUCTION 

Peripheral neuropathy is a common complication in 
diabetes mellitus (DM) that leads to serious functional 
disability; however, its evaluation has not been standar- 
dized (1-7). Electromyography with nerve conduction 
assessment is the neurophysiological test that confirms the 
distal symmetric polyneuropathy that is typical of diabetes, 
and it is the standard protocol in most centers that specialize 
in neuropathy (6-9). However, this test does not evaluate 
small fibers, is minimally invasive, requires expertise, and 
is time consuming (5,6). The clinical aspects, including 
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validated inventories, are valued in specialized diabetes 
centers; nevertheless, the screening for subclinical cases has 
low sensitivity, as demonstrated by the absence of symp- 
toms or signs of distal symmetrical polyneuropathy. 

Small fiber neuropathy may occur at any stage of 
diabetes, including during pre-diabetes. It is often painful, 
even in the absence of abnormalities, which are measured 
using conventional standard neurophysiological tests. The 
discrepancy between clinical presentation and test results 
could delay the correct diagnosis and appropriate treatment 
(8,10-12). Many diagnostic methods can assess small fibers; 
however, no single test can make this diagnosis. A battery of 
sensitive, reproducible and specific tests covering the 
somatic and autonomic systems is recommended (5,11,12). 
These tests include cardiovascular monitoring, sudomotor 
testing, pupillary responses, quantitative sensory tests, laser 
evoked potentials and thermography (12). 

Using such techniques has proven useful not only for 
diagnosis but also for guiding adequate therapy and 
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optimizing the follow-up. Infrared computerized thermo- 
graphy is a method of visualization, documentation and 
measurement of the infrared rays along the human body. 
According to Stefan-Boltzmann, this emission is propor- 
tional to the temperature of the skin and is directly related 
to blood flow in the cutaneous vessels (13,14,15). This 
method offers some advantages over others because it is 
completely non-contact, is painless and does not generate 
inputs beyond the equipment itself. 

In diabetes, vascular activity in the extremities, especially 
in the feet, may be indicative of autonomic neuropathy. 
Because vasomotor tone is regulated by fibers from the 
neurovegetative sympathetic system, its dysfunction could 
be associated with varying temperature patterns (16-19). 
Peripheral neurovegetative sympathetic nerve degeneration 
in advanced neuropathy damages the neurogenic control 
mechanisms that regulate capillary and arteriovenous (AV) 
shunt flow, which leads to an increase in the AV shunt in 
the feet of patients with diabetic neuropathy (20-23). These 
shunts are maintained in the constricted state by neurove- 
getative sympathetic tone. Losing this tone because of 
sympathetic neuropathy results in the shunt opening and 
the deviation of blood flow from the skin (20-23). 

Because of the great vulnerability that small fibers have to 
metabolic changes associated with hyperglycemia, injuries 
to these fibers may present earlier and with less marked 
changes than sensory-motor neuropathies (24,25). Although 
diabetic neuropathies are classified as diabetic microangio- 
pathic complications, it is now known that their pathophy- 
siological mechanism is multifactorial, and there is 
sufficient evidence that small-fiber polyneuropathy and 
even cardiac autonomic neuropathy may precede diabetes 
(i.e., during the pre-diabetic condition) (26-29). 

Cardiovascular reflex tests are standard for clinical 
autonomic evaluation because of their good sensitivity, 
specificity, and reproducibility; these tests are also non- 
invasive and well standardized. However, these tests 
require a well-trained evaluator and considerable evalua- 
tion time. Cardiovascular reflex tests can identify cardiac 
autonomic neuropathy even at the subclinical stage (29-36). 
This type of neuropathy can be subdivided into subclinical 
(predominantly functional and reversible) and clinical 
(structural neuronal changes are present). The subclinical 
subtype can only be diagnosed by tests and may occur as 
early as the time of the initial diabetes diagnosis. The clinical 
subtype is symptomatic and occurs during more advanced 
disease stages (35). 

The autonomic nerves are affected in various clinical 
diabetic neuropathy subtypes. In the most common type 
(typical polyneuropathy: symmetric, distal, and predomi- 
nantly sensory), there is a strong correlation between the 
progressive damage of the somatic and autonomic fibers. 
Indeed, 50% of diabetic patients with polyneuropathy have 
asymptomatic cardiac autonomic neuropathy, while 100% of 
patients with symptomatic cardiac autonomic neuropathy 
have polyneuropathy (3,37). The prevalence of cardiac 
autonomic neuropathy progressively increases in direct 
proportion to age, diabetes duration, and poor glucose 
control. 

The aim of this study was to evaluate plantar thermo- 
graphy sensibility and specificity in diabetic patients with 
polyneuropathy at the diagnosis using cardiac tests (heart 
rate variability) as a reference standard because the auto- 
nomic small fibers are first affected, and there is a close 
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relationship with cardiac and somatic autonomic nerve fibers 
in diabetes (18,23). 

MATERIALS AND METHODS 

This study was performed from March 2010 to August 
2011. The principles of the Declaration of Helsinki (38) were 
applied, and all of the patients provided informed written 
consent. The study was approved by the Ethics Committee 
of the Hospital de Clfnicas de Porto Alegre (HCPA), RS, 
Brazil, decision No. 09-446, January 2010. 

Subjects 

Seventy-nine individuals aged between 19 and 79 years 
(28 males) were evaluated and divided into three groups: 
controls (n = 37), pre-diabetics (n = 13), and type 2 diabetics 
(n = 29). The patients with pre-diabetes and diabetes were 
referred to the study from the Endocrinology and Pre- 
diabetes Out-clinic Unit of the Hospital de Clfnicas de Porto 
Alegre (Porto Alegre, Brazil). Diabetes and pre-diabetes 
were defined according to the American Diabetes 
Association (ADA) criteria (9). The control group consisted 
of volunteers that did not fit the criteria for type 2 diabetes 
or pre-diabetes. Smokers and subjects with other conditions 
that could potentially cause neuropathies, such as hypothyr- 
oidism and alcoholism, or a condition in which the thermal 
plantar images could mimic neuropathy, such as lumbosa- 
cral radiculopathy, were excluded. Ischemic heart disease or 
stroke patients were also excluded because of the heart rate 
variability tests were used. 

Data collection and instrumentation 

The data collection was conducted at a controlled 
laboratory temperature of 23 + 0.5°C and a relative humidity 
of 50 + 5%. All of the tests occurred on the same day and 
lasted 1.5 hours per participant. 

Clinical evaluation 

Clinical data, including age, gender, body mass index, and 
arterial systolic and diastolic pressure, were collected on the 
evaluation day. The clinical diagnosis of distal symmetric 
polyneuropathy was assessed using the Michigan inventory 
(39), a neuropathy score that consists of an inspection of foot 
deformities or ulcers and a brief neurological examination. We 
used the Achilles reflex test, a vibration sensitivity test on the 
hallux dorsum with a 128 Hz tuning-fork and a test for tactile 
sensibility using 10 g of nylon monofilament (Sory®, Bauru- 
SP-Brazil) on the plantar aspect of the hallux. The Michigan 
inventory was considered to be positive for neuropathy when 
four or more points (of a possible ten) were scored. 

Heart rate variability 

Both the ADA and the American Academy of Neurology 
(AAN) (40,41) recommend the following protocol: three 
Ewing tests (the deep breathing test, Valsalva maneuver, 
and orthostatic test) at the time of a type 2 diabetes 
diagnosis and five years after a type 1 diabetes diagnosis 
and repeated annually thereafter. These three tests per- 
formed together have good reproducibility and specificity 
above 91%; the deep breathing and orthostatic tests have 
93% sensitivity, and the Valsalva maneuver has 98% 
sensitivity. In addition, a spectral analysis was performed 
for heart rate variability, which is a result of sympathetic 
and parasympathetic balance at the sinus node. The method 
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consists of seven parameters that are evaluated together, 
including a three-band spectral analysis of heart rate 
variability and four Ewing tests (the Valsalva maneuver, 
orthostatic test, deep breathing test and orthostatic hypoten- 
sion test); this method has 98% sensitivity and 100% 
specificity (35). 

To evaluate the presence of cardiac autonomic neuropathy, 
heart rate variability tests were performed. These tests 
comprised three spectral indices (in the frequency domain) 
and four Ewing tests (35), including the Valsalva maneuver, 
orthostatic test, deep breathing test, and orthostatic hypoten- 
sion test. The electrocardiogram was recorded (particularly 
the QRS complex) using electrocardiography equipment 
(Neurosoft®, Ivanovo, Russia) and software that was created 
for heart rate variability analysis (Poly -Spectrum®). A ques- 
tionnaire concerning autonomic dysfunction symptoms was 
given following the protocol described by Rolim et al. (35). 

The presence of incipient cardiac autonomic neuropathy 
was defined as having two abnormal tests (98% specificity), 
and established neuropathy was diagnosed when three tests 
were abnormal (100% specificity) (35). In this study, we 
considered patients to be either positive or negative for 
cardiac autonomic neuropathy. 

Plantar thermography 

The plantar infrared images were recorded using three 
thermal infrared cameras (the PV320T Electrophysics, IRI 
4010 IRYSIS, and T400 Flir) with thermal sensitivities of 0.08 
and 0.01 °C, a full-scale spectral range between 7-12 and 8-14 
\i and a maximum error of 2%. The protocol for infrared 
images followed the recommendations of the American 
Academy of Thermology (42). Notably, caffeinated drinks or 
other vasoactive substances were suspended for at least 
eight hours prior to testing, and all prescribed medications 
were suspended 12 hours prior to testing. 

The participants were acclimated for 15 minutes in the 
examination room by lying on a stretcher with bare legs and 
feet and no surface contact. The plantar infrared image was 
recorded at baseline followed by the provocative maneuvers 
using the cold stress test. This test consists of immersing the 
feet, which are protected with thin plastic, for 60 seconds in 
water at 15 °C. After 10 minutes, a new plantar infrared image 
was recorded to evaluate the thermal recovery index. To 
calculate the recovery index, the average temperatures of 10 
regions of interest with similar dimensions were used: the 
hallux, 1st, 3rd, and 5th metatarsal heads and heel on both 
soles, as shown in Equation 1. 

\ TaB ( Lh " + L 1 mh + Limh + L5ml ' + Ll " + RJ <"+ R,ml <+ R i ml <+ R5ml >+ M"-' ) y 1 00l 

Y a YyCST Lha + Llmh + L3mli + L5mh + Lhe + Rha + Rlmh + lQmli + R5mh + Rhe 

TRI: thermal recovery Index; TaB: the average basal 
temperature of the 10 ROI ( C); TalO'CST: the average 
temperature 10 minutes after the cold stress test of the 10 
ROI ( C); L: left foot; R: right foot; ha: hallux; 1 mh: 1 st 
metatarsal head; 3 mh: 3 r metatarsal head; 5 mh: 5 th 
metatarsal head; he: heel. 

The interdigital anisothermal was assessed and was 
considered to be positive when thermal gradients (AT) 
>0.4°C existed between any of the toes 10 minutes after the 
cold stress test (41-49). Figure 1 shows two examples of the 
interdigital anisothermal test. 
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Figure 1 - (A) Plantar thermographic image in a diabetic patient, 
showing Interdigital Anisothermal (the white arrow shows the 
different colors in toes meaning AT > 0.4°C). (B) Plantar 
thermographic image in a control subject, with regular thermal 
distribution, without Interdigital Anisothermal (color gradient of 
toes considered normal, AT < 0.4°C). 



In this study, which is based on specific literature, an 
interdigital anisothermal result or a thermal recovery index 
of <90% or >100% after the cold stress test were considered 
(alone or in combination) to be positive for early neuropathy 
(small fibers and neurovegetative sympathetic fibers with 
vasomotor functions) (42-50). 

Electromyography 

The neuropathy screening protocol included the func- 
tional assessment of motor and sensory nerves of the four 
segments, as well as a myography using a needle electrode in 
suspected cases of axonal injury or root involvement (3,5-8). 
To record and analyze the data, a two channel electromyo- 
gram (Neurosoft®, Ivanovo, Russia) and dedicated software 
(NeuroMep®) were used. This method classifies the neuro- 
pathy into its various manifestations (i.e., mononeuropathies, 
multiple mononeuropathy and polyneuropathy). In this 
study, EMG was considered to be positive when it demon- 
strated distal symmetric polyneuropathy and was typical of 
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Table 1 - Patient characteristics. 



Variables* 



DM (n = 29) 



Pre-DM (n = 13) 



Control (n = 37) 



p-value 



Age (years) 
Gender (Female/Male) 
Time of diagnosis (years) 
Body mass index (kg/m 2 ) 
Systolic BP (mmHg) 
Diastolic BP (mmHg) 



55.9±9.4 b 
20/9 
5 (2-9) 
27.1 + 4.0 ab 
131.7±8.0 b 
85.9 + 7.9 b 



56.8±12.6 b 
10/3 
1 d-2) 
30.0±2.8) b 
130.5±6.2 b 
83.5 + 6.3 ab 



45.1 ±14.9 a 
21/16 

25.3(3.8) a 
121.7(9.0) a 
79.7(8.8) a 



0.001 
0.350 
<0.001 
<0.001 
<0.001 
0.012 



BP: Blood Pressure. * Mean +standard deviation or median (25 th percentile to 75 th percentile). 
5% significance level. 



' Same letters do not differ by Tukey's post hoc test at the 



diabetic neuropathy (5,7,51,52). The evaluator was blinded to 
the other test results. 

Statistical analysis 

The categorical variables were described by absolute and 
relative frequencies, and the continuous variables were 
described by the mean and standard deviation (SD) or the 
median and interquartile range. The cutoff point for infrared 
thermography, which was represented by the thermal 
recovery index in this study using heart rate variability as 
a reference standard, was obtained by creating a receiver 
operating characteristic (ROC) curve. The area under the 
curve was used to estimate the accuracy of the method. The 
Kappa coefficient was assessed to analyze the agreement 
between the thermographic methods versus heart rate 
variability. The association between the categorical variables 
was assessed using the chi-squared test or Fishers exact test. 
The means were compared using Student's t-test or a one- 
way analysis of variance (ANOVA) with Tukey's post-hoc 
test. In asymmetric cases, a Mann-Whitney U-test was 
applied. A significance level of 5% was adopted, and all of 
the statistical analyses were performed using SPSS 
(Statistical Package for the Social Sciences) version 18.0. 

RESULTS 

The analyses were separated by group because the groups 
were not homogeneous, with the exception of gender. The 
control subjects were younger {p- 0.001) and had a lower 
body mass index than the pre-diabetic subjects (p<0.001). 
The controls also had lower blood pressure values than the 
diabetic and pre-diabetic subjects (systolic, p<0.001, diasto- 
lic, p = 0.012), as shown in Table 1. 

For the three glycemic statuses, the pre-DM group 
presented with the highest thermal recovery index and 



Table 2 - Absolute and relative frequencies of neuropathy 
test results according to glycemic status. 

Test DM(n = 29) Pre-DM (n = 13) Control (n = 37) 

Abnormal Normal Abnormal Normal Abnormal Normal 

TRI (%) 22 (75.9) 7 (24.1) 10 (76.9) 3 (23.1) 14 (37.8) 23 (62.2) 

IDA(%) 20 (69.0) 9 (31.0) 10 (76.9) 3 (23.1) 13 (35.1) 24 (64.9) 

TRI+AID(%) 15 (51.7) 14 (48.3) 10 (76.9) 3 (23.1) 13 (35.1) 24 (64.9) 

HRV(%) 16 (55.2) 13 (44.8) 5 (38.5) 8 (61.5) 3 (8.1) 34 (91.9) 

Michigan (%) 2 (6.9) 27 (93.1) 0 (0) 13 (100) 0 (0.0) 37 (100) 

EMG (%) 16 (55.2) 13 (44.8) 2 (15.4) 11 (84.6) 0 (0.0) 37(100) 

TRI: thermal recovery index; IDA: interdigital anisothermal, HRV: heart 
rate variability, (Michigan inventory); EMG: electromyography. 



had abnormal interdigital anisothermal results. The DM 
group presented with the highest number of abnormal 
results from the neuropathy diagnosis tests (heart rate 
variability, Michigan inventory, and EMG), as reported in 
Table 2. Although we observed a greater number of positive 
results for neuropathy diagnosis, only 6.9% (n = 2) of the 
patients from the DM group presented with a positive 
Michigan inventory test, thus identifying clinical neuro- 
pathy by the loss of protective foot sensation. 

Cardiac autonomic neuropathy (diagnosed using the 
heart rate variability method) was the reference used to 
estimate the sensitivity and specificity of the plantar 
infrared thermography. The interdigital anisothermal test 
and the thermal recovery index, alone or in combination, 
were used to identify the subjects with cardiac autonomic 
neuropathy in the three groups (Table 3). In the diabetic 
subjects, the interdigital anisothermal test alone performed 
better than the thermal recovery index, and both had better 
sensitivity (81.5%) and specificity (46.2%). For the pre- 
diabetic patients, all three tests performed equally well 
(sensitivity 80.0%, specificity 25.0%). None of the control 
subjects had abnormal results for the interdigital anisother- 
mal test or the thermal recovery index, which would have 
precluded the sensitivity estimation in this sample of 
subjects. However, the specificity was higher in the control 
subjects (70.6%). 

Because the interdigital anisothermal test performed best, 
the clinical data regarding the results of this test are 
depicted in Table 4. In the pre-diabetic patients, the clinical 
data were significantly different between the subjects who 
presented or did not present with interdigital anisothermal 
results. The subjects in the pre-DM group who presented 



Table 3 - Comparison of the methods. 

Comparisons % abnormal Sensibility (%) Specificity (%) 



DM Group (n = 29) 



HRVxTRI 


55.2x75.9 


75.0 


23.1 


HRVxIDA 


55.2x69.0 


81.3 


46.2 


HRVxboth (TRI+IDA) 


55.2x51.7 


56.3 


53.8 


Pre-DM Group (n = 13) 








HRVxTRI 


38.5x76.9 


80.0 


25.0 


HRVxIDA 


38.5x76.9 


80.0 


25.0 


HRVxboth (TRI+IDA) 


38.5x76.9 


80.0 


25.0 


Control Group (n = 37) 








HRVxTRI 


8.1 x37.8 




58.8 


HRVxIDA 


8.1 x35.1 




61.8 


HRVxboth (TRI+IDA) 


8.1 x27.0 




70.6 



HRV: heart rate variability; TRI: thermal recovery index, IDA: interdigital 
anisothermal. 
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Table 4 - Interdigital anisothermal clinical characteristics. 



Variables* 




DM (n = 29) 




Pre-DM (n = 13) 




Control (n = 37) 






Abnormal 


Normal 


p-value 


Abnormal 


Normal 


p-value 


Abnormal 


Normal 


p-value 


Age (years) 


56.1 ±8.8 


55.4±11.1 


0.876 


60.7 ±10.9 


43.7±9.1 


0.032 


46.4±16.4 


44.3 ±14.3 


0.695 


Gender (Female/ 


14/6 


6/3 


1.000 


9/1 


'A 


0.108 


8/5 


12/12 


0.744 


Male) 




















Time of diagnosis 


5.5 (2.3-9.5) 


2 (1-10.5) 


0.216 


1.5 (0.9-2.0) 


1 (1-D 


0.469 








(years) 




















BMI (kg/m 2 ) 


27.2 + 4.4 


26.9±3.3 


0.825 


29.7±3.0 


31.1 ±2.2 


0.467 


26.2 ±4.6 


24.8±3.2 


0.290 


SBP (mmHg) 


133.5±5.9 


127.8 ±10.9 


0.170 


129.7±6.7 


133.3±2.9 


0.393 


121.2±10.8 


121.9±8.1 


0.799 


DBP (mmHg) 


87.0±5.5 


83.3± 11.7 


0.393 


83.5 ±11. 7 


83.3±7.6 


0.970 


80.0 ±11.7 


79.6±7.1 


0.893 


*Mean (SD) or median (25 th percentile to 75 th percentile). BMI: body mass index; SBP: systolic blood 


pressure; 


DBP: diastolic blood pressure. 


Level of 



significance = 5%. 



with an abnormal plantar thermography result were older 
than those who presented with a normal interdigital 
thermographic pattern. 

DISCUSSION 

The presence of interdigital anisothermal using a simpli- 
fied plantar thermographic study appears to be the most 
appropriate diagnostic test for the early diagnosis of 
neuropathy in the patients with DM and pre-DM. Because 
of its sensitivity and specificity, in addition to its conve- 
nience and timing, interdigital anisothermal assessed by 
plantar thermography seems to be most suitable for diabetic 
neuropathy screening programs. Its application is simple 
using the proposed protocol, and it is cost-effective because 
it requires no additional supplies. 

It is important to identify diabetic neuropathy early to 
prevent secondary complications, such as neuropathic pain 
and diabetic foot (10,11,52,53). Unfortunately, because of 
economic concerns and the lack of technological resources, a 
consensus of international associations dedicated to treating 
and preventing diabetes has advocated only using clinical 
examination as a population screening method (56-59). In 
this study, the clinical application of a validated clinical 
inventory (Michigan inventory) did not identified neuro- 
pathy even in cases with positive EMG, which would 
certainly include late cases of neuropathy. Plantar thermo- 
graphy is a new, non-invasive method that can be included 
in neuropathy screening programs, thereby increasing the 
sensitivity for proper diagnosis. 

The current method for the early diagnosis of diabetic 
neuropathy is a battery of tests, including methods capable 
of assessing thick and thin nerve fibers, particularly the 
autonomic nerve fibers (5,12). The two types of thermo- 
graphic tests (thermal recovery index after the cold stress 
test and the interdigital anisothermal test) demonstrated 
differences in patterns between the three groups. These 
patterns reflect the sympathetic nerve fibers with vasomotor 
tone function and its effect on the skin temperature of the 
feet (15,17,18). 

Another benefit provided by plantar thermography for 
chronic patients, such as those with DM, is the possibility of 
following the evolutionary aspects of the disease because 
thermography allows for functional imaging (60). Because 
of the strict control of risk factors, such as glycemic control, 
cutaneous vasomotor functionality of the feet can be 
monitored periodically using this method. 



In this study, plantar thermography proved useful in the 
early diagnosis of diabetic neuropathy, particularly the 
small and autonomic fibers that are commonly associated 
with a sub-clinical condition. Future research studying the 
thermal recovery index for diabetes duration would be 
useful, particularly when comparing the density of nerve 
fibers using skin biopsy and Quantitative Sensitive Test 
(QST) to demonstrate a low sensitivity to heat. In addition, it 
would be interesting to compare the potential amplitude of 
the sural nerve, when present, with the thermographics 
variable studied here (interdigital anisothermal and thermal 
recovery index), as demonstrated by Shun et al. (61). 

One limitation of this study is that it is a cross-sectional 
study that involves a sensitive test. Therefore, follow-up will 
be needed to assess the development of more advanced 
forms of neuropathy and diabetic complications. 
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